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SUMMARY 


After  a  discussion  of  the  basic  shortcoming  of  the  mathematical 
expression  generally  used  to  compute  atmospheric  refraction  in  con¬ 
nection  with  satellite  observations,  an  approach  is  presented  for 
computing  corrections  for  refraction  for  both  the  geodetic  satellite 
triangulation  and  the  topographic  photogramme trie  satellite  method. 

For  target  points  outside  the  effective  atmosphere,  the  refrac¬ 
tion  is  obtained  as  astronomical  refraction  minus  a  correction  angle. 
This  correction  angle  is  a  function  of  the  corresponding  astronomical 
refraction  and  at  the  same  time,  constitutes  the  amount  of  refraction 
encountered  if  a  photogrammetric  camera  i6  placed  outside  the  effec¬ 
tive  atmosphere. 

This  correction  angle  is  insensitive  to  changes  in  astronomical 
refraction.  It  is,  therefore,  concluded  that  the  determination  of 
refraction  from  a  ground  station  to  a  satellite  is  mainly  affected 
by  the  error  made  in  determining  astronomical  refraction,  and  direc¬ 
tions  observed  at  a  satellite  are  affected  only  insignificantly  by 
refraction  anomalies.  In  addition,  approximation  formulas  are  given 
for  use  in  connection  with  moderately  sized  zenith  angles. 


iv 


THE  INFLUENCE  OF  ATMOSPHERIC  REFRACTION 


ON  DIRECTIONS  MEASURED  TO  AND  FROM  A  SATELLITE 


I.  INTRODUCTION 

The  geodetic -photogrammetric  community  is  considering  a  twofold 
purpose  for  satellites.  For  supporting  classic  geodetic  triangula¬ 
tion  schemes,  the  satellite  is  used  as  an  auxiliary  target  point  and 
the  direction  to  the  satellite  is  measured  from  the  ground.  For 
topographic  mapping,  the  satellite  serves  as  the  carrier  for  a  pre¬ 
cision  photogrammetric  camera  system — measuring  the  direction  from 
the  satellite  to  the  ground.  For  both  objectives,  the  accuracy  re¬ 
quirements  are  such  that  the  problem  of  atmospheric  refraction  must 
be  considered. 


II.  GENERAL  PRINCIPLES  OF  THE  REFRACTION  GEOMETRY 

Solutions  for  the  numerical  evaluation  of  the  problem  have  been 
published  in  the  past  [Brown,  1957;  Case,  1962;  Holland,  1961; 
Jones,  I96I;  and  Schmid,  19591*  These  solutions  and  this  one  are 
based  on  the  classic  refraction  geometry  as  shown  on  page  2.  A  rep¬ 
resentative  expression  for  these  solutions  is: 


tan  za 


sin  Cs 


cos  Cs  - 


£a 

rs 


(1) 


where 


Cs  =  (z)a  +  ^  -  Ar  and  Ar  =  sin-1 

rs 

k  =  nara  sin(z)a  r*,  denotes  the  astronomical 

refraction 

Correctly,  the  practical  value  of  an  expression  of  the  form  of  form¬ 
ula  (l)  has  been  questioned.  The  strength  of  such  a  solution  has 
been  doubted  because  of  the  usually  small  central  angle  C8  (e.  g. 
[Case,  1962]).  More  correctly,  the  weakness  of  these  solutions  is 

rs 

caused  by  the  ratio:  - - - 

rs  "  ra 

A  differentiation  of  formula  (l)  with  respect  to  Cs  gives: 


1 


2 


(2) 


— 3 — 

cos  za 


Az; 


rs 

r.aJ?.°ldJ 


-  1 


a  "  rs  cos  Cs  *a 

ra  +  rs  G0S  cs 


AG£ 


-  2 


if  Cs  is  a  small  angle  cos  Cs  «  1.  Thus,  one  obtains: 


or 


or 


(3) 


(4) 


(5) 


where  Hs  is  the  height  of  the  orbiting  satellite. 

It  follows  from  formula  (5)  that  formula  (l)  becomes  the  more 
useful  the  greater  the  height  of  the  orbiting  satellite.  Correspond¬ 
ing  error  coefficients  are  given  in  Table  I. 

Table  I.  Numerical  Evaluation  of  Formula  (5) 
for  Aza  for  Various  Heights  of  Satellite 


Zenith 

Angle 

(0) 

Hs  =  100  km 

(sec/sec) 

Hs  =  300  km 

(sec/sec) 

Hs=  1,000  km 

(sec/sec) 

Hs  =  ra 

(sec/sec) 

Hs  =  10  ra 

(sec/sec) 

0 

64.70 

21.00 

7.00 

2.00 

1.10 

15 

60.37 

19.59 

6.53 

1.87 

1.03 

30 

48.52 

15.75 

5.25 

1.50 

0.82 

45 

32.35 

10.50 

3.50 

1.00 

0.55 

60 

16.18 

5.25 

1.75 

0.50 

0.28 

75 

4-33 

1.4l 

0.47 

0.13 

0.07 

However,  both  the  unfavorable  error  propagation  and  the  neces¬ 
sity  of  computing  with  the  small  central  angle  Cs  can  be  avoided 
without  losing  the  strength  of  the  geometry  inherent  in  a  spherical 
earth  solution. 
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(6) 

(7) 


From  the  classic  refraction  geometry  (page  2): 
za  =  (z)a  +  ra  =  (z)a  +  -  a 

and  zs  =  (z)s  +  rs  =  (z)s  +  a 

From  formulas  (6)  and  (7),  it  follows  directly  that: 

ra  +  rs  =  r^  (8) 

A.  Directions  Measured  from  Ground  to  Satellite  (Spherical  Earth). 
Again  directly  from  the  classic  refraction  geometry: 


k  =  nr  sin  (z)  =  nara  sin  (z)a 

(9) 

Ak  =  ra  sin  Cx  =  ra  sin  [ (z)a  +  r,J 

(10) 

and 

k*  =  ra  {  na  sin(z)a  -  sin  [  (z)a  +  r«,]  1 

(n) 

,  -1  k*  k* 

a  =  sin  -3-  m  -3-  p 

a  a 

(12) 

where  d  is  the  distance  between  ground  and  satellite.  (All  linear 
parameters  are  in  meters.) 

A  classic  expression  [Jordan,  19391  for  the  relation  between 
astronomical  refraction  (rM )  and  index  of  refraction  (n)  is : 

roo  =  (n-l)  (l-s)  tan(z)  -  |^s(n-l)  -  ( n J  tan3  (z)  (13) 

where  s  =  F  —  =  0.001255  (lU) 

g  *•(!) 

P  =  weight  of  mercury  column 
O.76  m  length  where  the 
density  of  mercury  =  13.596 

(S) =  density  of  air  =  0.001293 

r  =  6,370,000  m 

which  may  be  written  as: 


k 


RT, 


s  = 


0.001255 


(15) 


where  R  =  gas  constant  =  2 9. 27*+ 5 
Tq  =  273.16  K° 

Formula  (13)  can  be  rearranged  to  read: 

r=c  =  (na-l)  tan ( z ) a  £l  -  s  -  tan2  (z)a(s 
A  first  approximation  for  (na-l)  is  obtained  with 


■¥5 


(16) 


na"1  tan(z)£ 


(17) 


Substituting  formula  (17)  into  the  parenthetic  terms  of  formula  (16), 
one  obtains: 


Qa  ^  +  I  tan(z). 


where 


1  =  1 


roo  tan(z). 


Uos2(z)£ 


(18) 


(19) 


Substituting  formula  (l8)  into  formula  (ll) 

k*  =  ra  j^sin(z)a(l  -  cos  r„)  +  cos(z)a  -  sin  r^j  (20) 

Because  r,*,  is  a  small  angle,  sin  rM  •*  r*,.  Therefore,  formula  (20) 
can  be  written  as: 

k*  =  ra  Qin(z)a(l  -  cos  r^)  +  cos(z)a  -  l^J  (21) 

Because  the  parenthetic  term  in  formula  (19)  is  small,  one  can  write: 
I  =  1  -  x  and  consequently  ^  -  3^  =l+x-l=x 


where 


_  s  _  rac  tan(z)a 
cos^(z)a  2 


(22) 
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With  the  corresponding  substitution,  formula  (21)  reduces  to: 
k*  »  ra  j(sin(Z)a(l  -  cos  rj  +  ] 


which  can  be  arranged  into  the  form: 


k*  =  ra  sin  (z) 


a 


E‘- 


cos  O  - 


0 


rar9cs 

cos  (z). 


r^  being  small,  cos  *.  1  — and,  therefore, 


[ 


(1  -  cos  r  ) 


Thus,  one  obtains  finally: 


-#]■ 


(23) 


(2k) 


(25) 


r.r  S 

k*  = 

cos(z j 


(26) 


sind  with  formula  (15 ) 

*■  -  BTc(l  *  %)  5^  <»a  8eing  <  8,000  m)  (27) 

With  formula  (12),  the  satellite  refraction  a  is  obtained  as 

(28) 


a  = 


*ar~  *  3 


=  d  cos(z). 


H„ 


ra  *  s  *  roo 

"F  Hs  tan2(z)a~j 

’  12.5  ’  lO^J 


Z\-  Ha 


Because  of  the  insensitivity  of  this  expression  with  respect  to  d, 
for  moderately  sized  zenith  angles  it  is  possible  to  set 


Hs  -  Hg, 

d  -  ,  y-  (flat  earth  approximation) 

COS l Z  } 

v  'a 


(29) 


Consequently:* 


rn  s 


vS 


*■"  -  RT°  TSTHij 


(30) 


*  For  limitations  of  formula  (30),  compare  results  in  Table  IV. 
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In  order  to  study  the  propagation  of  errors  of  ra>  ,  and  (z)a  with 
respect  to  k*,  formula  (26)  is  differentiated: 


5k*  =  _ 

~  cos 


So  8 

(z)a  a 


k* 


<Jr« 


5k*  s  tan  (z)a  5(z)a  ^  <?(z)a 

“  tt  ®  tan yZ Jq  ti 


5k* 


cos 


ra  s 
cos(z) 


TzTi 

T 


k* 


and  accordingly  with  formulas  (28)  and  (30) 


=  JL.  5r  and 

***.  ra  a 


da 


da 

5^51 


=  a 


tan(z). 


5(HS  -  Ha) 

<5(z)a 


He 


— -  5(H  -  H  ) 

Ho  '  s  a' 


(31) 

(32) 

(33) 

m 

(35) 


4^’  ») 

Mean  astronomical  refraction  can  he  computed  according  to  de  Ball 
[  1906]. 

r~m"  *6  tan(z)a  "  Ai  tan3(z)a+  *2  tan5(z)a-  A3  tan7(z)a  +  . . .  (37) 


A  similar  formula,  useful  for  extremely  large  zenith  distances,  was 
developed  by  Garfinhel  [19*^] • 

de  Ball's  formula  was  used  in  Table  II  which  gives  in  column 
(5)  mean  astronomical  refraction  room  for  various  zenith  distances. 

Astronomical  refraction  rM  is  then  obtained  as: 

•  W  (38) 

where  W  is  a  meteorological  correction  factor.  For  comparison, 
column  (6)  in  Table  II  gives  astronomical  refraction  for  t  =  10°  C. 

w  _  *a  (1  +  foo) 

P0  1  +  0ta 

PQ  =  76O  mm/Hg 
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t0  =  0°  c 
0  =  0.003665 


Table  II.  Mean  Astronomical  Refraction 


~rrr~ 

Zenith 

Angle 

(*« 

—J2) 

A 0  tan(z)a 

II 

- T3T - 

-A{  tan3(z)a 

It 

“ 1*1 - 

A^>  tan-3(z)a 

II 

(5)  '■  ' 
r 

80  m 

n 

- (5) 

for  ta  =  10°  C 

It 

0 

0 

0 

0 

0 

0 

15 

16.104 

-0.001 

0 

16.10 

15-53 

30 

34.669 

-0.013 

0 

34.69 

33-46 

45 

60.101 

-0.067 

0 

60.03 

57.90 

60 

104.097 

-0.346 

0.003 

103.75 

100.08 

75 

224.298 

-3*461 

0.153 

220.99 

213.27 

Table  III.  Numerical  Evaluation 
Formulas  (20),  (26),  (31),  (32),  (33), 

Computed  with  r  =  r  and  H„  = 
“  “m 

of 

and  (4l) 

0 

(1) 

Zenith 

Angle 

(z)a 

(2) 

Formula 

(20) 

(3) 

k*(m) 

Formula 

(26) 

(4) 

Formula  (4o) 
or  (4l) 

(5) 

<?k* 

tea. 

m/m 

(6) 

dk* 

<?(z)a 

m/sec 

(7) 

<?k* 

dr 

m/sec 

0 

0 

0 

0 

0 

0 

0 

15 

O.67 

0.65 

O.65 

0.0000001 

0.000001 

0.040 

30 

1.55 

1.55 

1-55 

0.0000002 

0.000004 

0.045 

45 

3.26 

3.29 

3.29 

0.0000005 

0.000016 

0.055 

60 

8.07 

8.04 

8.06 

0.0000013 

0.000068 

0.077 

75 

33.54 

33.07 

33.56 

0.0000052 

0.000598 

O.I58 

From  column  (7)  of  Table  III,  it  is  seen  that  the  computation  of  k* 
is  insensitive  to  an  error  in  r*  .  Table  II  shows  that  the  first- 
order  term  of  formula  (37)  gives  a  good  approximation  for  rM  .  Con¬ 
sequently,  formula  (26)  can  be  written  as:  m 


k* 


tan  (z)a 
cos  (z)a 


W 


or  with  formula  (27)  as: 


(39) 


8 


with 


.  tan  (z)a 
k*-RIoAi  cos  [»)* 


.  w 


(ho) 


r  =  6,370,000 

rt0  =  7996.622420 
=  O.OOO29137566 


k* 


2.330 


tan(z)£ 
COB  (z  ), 


W 


(41) 


The  adequacy  of  formula  (4l)  can  be  Judged  from  Table  III  by  compar¬ 
ing  columns  (2)  and  (4). 

Again  with  formula  (12),  one  obtains  from  formula  (4l) : 
a" 

and  in  accordance 

a  " 

For  limitations  of  formula  (43),  compare  corresponding  results  in 
Table  IV. 


2.3300  tan  (z)a 
d  cos  (z)a 


p"  •  W 


(42) 


with  formula  (29) : 


-  1^20  tan  (z) 


Hs  -  Ha 


W 


(43) 


t v  le  IV  shows  the  computation  of  a  for  various  zenith  angles 
and  heights  of  the  orbiting  satellite  with  formulas  (28),  (30),  (42), 
and  (43)  for  Ha  =  0.  Table  V  shows  the  corresponding  error  coeffi¬ 
cients  for  formulas  (34),  (35);  and  (36)  using  a  computed  with  form¬ 
ula  (42). 
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Table  IV.  Numerical  Values  for  a  for  Various  Zenith  Angles  and  Heights  of 
Satellite  Computed  with  Formulas  (28),  (30),  (42),  and  (43) 


H 

II 

a. 

W 


o 

s 


CQ 

W 


UJ 

1 


S'1 


J 


o 

o 


at 

W 


(2  OJ 
•H 

o 

T)  CO| 
4)  •>— 

I- 

a  co 

O  CM 

O 


a 


iS  OJ 
•H 


2 


a  co 

O  CVJ 
o 


rrt\ 

co  -4 

1 


O  OJ 

pt)  -4 

$ 


9* 


& 


O  co 

O  OJ 


5  0) 

•H  r-i 

a  5P 

as 


coco  co  co  o\ 

H  CVJ  -4  CO  t- 

•  •  •  •  • 

O  O  O  O  O  H 


COCO  H  CD  CO 
H  OJ  UNCKaD 

O  O  O  O  O  OJ 

coco  cocot— 
H  oi  -4  co  t- 

•  •  •  •  • 

O  O  O  O  O  H 


coco  i— I  co  t— 
H  OJ  unOncO 

•  •  •  •  • 

O  O  U  O  O  C\J 


ro 


as 


O  O  O  H 


CO 

-4 


toj 

On  no 


C^ao 

r—  cFn 

•  • 

OJ  LT\ 


ir\co 

ON-S 


O  O  O  H  OJ  t- 


on 


as 


O  O  O  H 


NO  On 

C— ao 

OJ  un 


co 

-4 


(04 

ON  NO 


UN  t— 

ON -4 


O  O  O  H  OJ  t- 


ON  t'-  H 
OJ  t“-  CO 


OJ  C 

coc 


OH01403N 


ONCO  -4 
OJ  t~-CO 


r-l  O 
UN  C~- 


O  H  OJ  -4  00  ON 

rH 


On  O 
OJ  t— co 


ON  t- 
OJ  NO 


O  H  OJ  -4  CO  t— 

rH 


ON  CO  .4 
OJ  t“C0 


o 

UN  NO 


O  H  OJ  -4  CO  ON 

H 


O  lOO  UNO 
H  (04  'Z> 


aj 

Ph 

o 

(— I 


at 

W 


M, 

W 


CO| 

CO  -4 


aJ- 

I 


o  OJ 
pc*  -d- 

•rH 

-d  co| 
<u  — - 

% 

O  co 

O  OJ 


(fl / 

IS1 


f§  Of- 

•rl 

>  O' 
-c)  co| 

d) 

£ 

&C0 
O  OJ 

o  — « 


43 

-p  aj 
•H  i— I 

a  «)  o 

at 


o  o  o  o  o  o 


OJ  UN  O  UNOI 
OOrlOJO 
O  O  O  O  H 

*  •  I  *  t 

o  o  o  o  o  o 


QJ  -4  CO  coco 

O  O  O  rl  OJ 
O  O  O  O  O 

•  •  •  •  • 

o  o  o  o  o  o 

OJ  LTN  O  LTN  CVJ 

88388 

•  •  •  •  • 

o  o  o  o  o  o 


01-4  co  coco 
o  o  o  r-j  oj 

o  o  o  o  o  o 

OJ  UN  on  o  CO 

OOOOIt- 

•  •  •  •  • 

o  o  o  o  o  o 


OJ  -4  CO  COCO 

O  O  O  rH  OJ 

•  •  •  •  • 

o  o  o  o  o  o 


OJ  UN  ON  O  co 
OOOOIP 

•  •  •  •  • 

o  o  o  o  o  o 


O  UN  o  UN  O  UN 
H  CO  -4  NO  C 
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Table  V.  Error  Coefficients 


or 

(2) 

— m — 

—(V - 

- (5) - 

- (5) - 

“771“ 

(Bl 

Zenith  da 

Angle  dra 

da 

<?HB 

& 

0 

(sec/m) 

(sec/m) 

(sec/sec) 

Ha  -  100  km 

3S  -  100  km 

Hb-  300  km 

Hr  - 1000  km 

Hfi-ra 

H«-  10ra 

0 

0 

0 

0 

0 

0 

0 

0 

15 

0.0000002 

0.0000129 

0.080 

0.027 

0.0081 

0.0012 

0.00012 

30 

0.0000004 

0.0000278 

0.080 

0.027 

0.0081 

0.0014 

0.00014 

45 

0.0000008 

0.000048'i 

0.081 

0.027 

O.OO85 

0.0015 

0.00017 

60 

0.0000013 

O.OOOO85I 

0.082 

0.028 

0.0094 

0.0019 

0.00024 

75 

0.0000031 

0.0001970 

O.O89 

0.034 

0.0130 

0.0033 

0.00046 

B. 

Directions 

Measured 

from  Satellite  to 

Ground  (Spherical  Earth). 

From  the  illustration  (page  2): 

(z)a  =  180  -  [(z)s  -  Cs  +  rj  (44) 


(ra  +  ®r) 


rs  sin  (z)s 
sin[(z)s  -  Cs] 


(45) 


where  the  refractive  height 


sin  t(z)s  -  Cs] 

Consequently, 

rs  sin(z)s  -  k* 
&  sin[  (z)s  -  CS1 


(46) 


(47) 


Because  of  the  insensitivity  of  the  computation  of  k*  to  a  change  in 
ra,  it  is  adequate  to  write: 


rg  sin  (z)s 
sin  t(z)s  -  CS1 


(48) 
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Substituting  formula  (44),  neglecting  ,  and  substituting  formula 
(48)  into  formula  (26): 


k*  = 


-2rs^  *  s  sin  (z)e 


sin  2  [(z)s  -  Cs] 

The  differentiation  of  formula  (49)  gives : 

5k*  k*  , 

=  —  5r 
rs 

5k*  k*,,  „  „r/\  _n. 

=  -pr  Utn  (z)s  -  2  ctn  2[(z)s  -  Cj  I 

3s;  =  7  2  Ctn  2  t(z)s  -  CS1  -jfr 
5k*  k*  5r» 

rM  p 


(49) 


(50) 

(51) 

(52) 

(53) 


In  practice,  instead  of  using  formula  (49),  it  will  be  from  the 
computational  standpoint,  advantageous  to  compute  with  formula  (44) 
neglecting  r,,,,. 


(z)a  =  180  -  [(z)s  -  CS1 
k*  is  then  computed  with  formula  (26)  as: 


(54) 


k*  = 


raroo  s 


cos  [ (z ) s  -  Cs] 
or  with  formula  (4l)  as: 


(55) 


tan  [(z)n  -  CJ 

k*  =  2.3300  - -7-fS - -f—  W 

cos  [(z)s  -  Cs] 

Using  formula  (56),  one  obtains  with  formula  (42) 


(56) 


,»  =  2.3300  tan  [(z)H 
d  cos  [(z)s 


^■p”  •  W 


-  Cgl 


(57) 


or  in  analogy  to  formula  (43) 
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(58) 


a"  =  tan  [(z)8  -  C81  p"  •  W 

For  limitations  of  formula  (58),  compare  the  corresponding  re¬ 
sult  in  Table  VI. 

Formulas  (49),  (56),  (57),  and  (58)  are  numerically  evaluated 
in  Table  VI.  The  result  is,  as  it  obviously  has  to  be,  in  agreement 
with  the  corresponding  results  presented  in  Tables  III  and  IV. 


Table  VI.  Numerical  Evaluation  of 
Formulas  (49),  (56),  (57),  and  (58) 


— m — 

- (3) - 

hoh 

- (57 

HI 

- (77 

Zenith  Angle 

k*  (m) 

a " 

Hr  (m) 

(z)a 

t(z)s-cs1 

Formula 

Formula 

Formula 

Formula 

Formula 

(49) 

(56) 

(57) 

(58) 

(46) 

Hs  =  100  km 

0 

180 

0 

0 

0 

0 

0 

15 

165 

O.65 

O.65 

l.?9 

1.29 

2.49 

30 

150 

1.55 

1*55 

2.78 

2.77 

3-10 

45 

135 

3.29 

3*29 

4.84 

4.81 

4.65 

60 

120 

8.03 

8.06 

8.52 

8.33 

9.27 

75 

105 

33*08 

33*59 

19*71 

17.94 

34.24 

c. 

Formulas  Necessary  for  Computing  Directions 

to  and 

from  a 

Satellite . 


Corrected  for  refraction,  the  formulas  necessary  for  computing 
directions  to  and  from  a  satellite  are: 


za  =  (z)a  +  5.  “  17 
zs  =  (z)s  +  a 


(6) 

(12) 


a  =  2 « 33  tan 
d  cos 


(z)a  p"  .  w  =  tanKz)s-c!&1  # 

vz)a  P  d  cos[(z)a-C«] 


s^sJ 


W 


(42) 

and 

(57) 


or 


Hi"  tan(z^a  ‘  P"  ■  W  =  tan t(z)s-Cs]  •  p"  .  W 

s  a  s 


(43) 

and 

(58) 
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(z)a  and  (z)6  are  the  measured  directions  at  the  ground  or  at  the 
satellite,  respectively.  Because  a  is  a  small  angle,  only  a  few 
digits  must  he  carried  in  the  corresponding  computations. 

can  he  computed  with  formula  (37)  or  according  to  Garfinkel 
[1944],  W  is  obtained  from  formula  (38). 


III.  DISCUSSION 

Formula  (6)  shows  that  a  direction  measured  on  the  ground  to  a 
satellite  must  be  corrected  for  astronomical  refraction  and  a  cer¬ 
tain  parallactic  correction  a  which  is  functionally  related  to  as¬ 
tronomical  refraction.  A  direction  measured  at  a  satellite  needs 
only  to  be  corrected  for  the  same  parallactic  correction  a  (formula 

(7)). 


The  computation  of  astronomical  refraction  is,  as  Table  I  dem¬ 
onstrates,  sufficiently  proportional  to  tan  (z)  (flat  earth  approxi¬ 
mation)  for  moderately  sized  zenith  distances,  and  for  larger  zenith 
distances,  higher  order  terms  must  be  carried  in  order  to  satisfy 
the  accuracy  requirements  for  geodetic  satellite  triangulation  pur¬ 
poses  .  In  practice,  it  will  be  necessary  to  either  determine  rM 
independently  by  elevation  angle  measurements  to  stars  or  by  in¬ 
directly  considering  the  astronomical  refraction  by  interpolating 
the  photographed  satellite  image  into  the  background  of  the  surround¬ 
ing  star  imagery. 

The  computation  of  k*  and,  correspondingly,  of  the  angle  a  is, 
as  shown  in  Tables  III  and  V,  extremely  insensitive  to  a  change  in 
r,  similarly  insensitive  to  changes  in  (z)a  of  t(z)s  -  Cs]  respec¬ 
tively,  and  fairly  insensitive  to  changes  in  r^,. 

Consequently,  satellite  photogrammetry  will  not  encounter  any 
difficulties  in  computing  the  corresponding  refraction  corrections. 
Because  a  is  a  small  angle  (data  in  Table  IV  or  [1959  ])>  it  is  feas¬ 
ible  to  first  compute  the  photogrammetric  triangulation  without  any 
refraction  correction.  Such  a  result  certainly  will  assure  more 
than  the  necessary  fidelity  in  the  geometric  relation  between  ground 
and  satellite,  for  computing  the  corresponding  refraction  corrections 
in  the  last  iteration  cycle. 

In  order  to  give  some  idea  of  the  error  caused  by  neglecting 
refraction  for  satellite  photogrammetry,  the  formula  (46)  for  re¬ 
fractive  height  Hr  (page  2)  is  given  and  corresponding  values  are 
presented  in  Table  VI,  column  7»  In  practice,  t(z)s  -  Cs]  will  be 
>  120°  resulting  in  an  Hr  <  9m*  Assuming  an  error  in  elevation  of 
about  1/20,000  of  flying  height  as  a  practical  limit  in  photogram¬ 
metric  triangulation,  it  becomes  evident  that  for  satellites  orbiting 
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at  200  to  300  km  height,  it  is  hardly  necessary  to  consider  refrac¬ 
tion  at  all. 


Concluding,  it  should  be  mentioned  that  the  results  obtained  by 
Jones  [1961]  have  given  rise  to  wrong  conclusions.  His  error  coeffi¬ 
cients  have  no  physical  significance,  but  are  essentially  the  values 
obtained  from  formula  (5)  in  Table  I. 

Holland  [1961]  is  considering  the  "Effects  of  Atmospheric  Re¬ 
fraction  on  Angles  Measured  from  a  Satellite."  Qualitatively,  he 
obtains  a  correct  conclusion,  but  his  approach  lacks  geometrical 
strength  and,  therefore,  his  formulas  and  numerical,  results  are  not 
adequate  for  computing  either  satellite  refraction  or  corresponding 
error  propagation  coefficients. 

As  seen  from  formula  (ll),  k*  is  the  difference  between  two 
large  numbers.  It  is,  therefore,  necessary  to  assure  geometrical 
fidelity  when  the  relation  between  the  index  of  refraction  and  the 
corresponding  astronomical  refraction  is  introduced.  Therefore,  for 
this  phase  of  the  problem  the  use  of  flat  earth  expressions  is  not 
adequate.  Holland  obtains  values  for  the  satellite  refraction  and 
for  the  corresponding  error  propagation  coefficients  which  are  ap¬ 
proximately  one  to  two  magnitudes  too  small. 


Case  [1962]  correctly  points  out  that  the  basic  formula  (l)  is 
not  suited  for  numerical  work.  However,  he  attributes  the  difficul¬ 
ties  to  the  small  center  angle  C8  and,  consequently,  supports  Jones' 
misleading  conclusion  concerning  the  propagation  of  an  error  in  as¬ 
tronomical  refraction.  By  using  Holland's  extended  formula  for  as¬ 
tronomical  refraction,  Case  obtains  realistic  values  for  a.  A  com¬ 
parison  of  his  results  in  Table  I,  column  (5),  with  Table  IV  of  this 
paper  shows  that  Case's  values  are  slightly  too  large  for  large 
zenith  angles,  as  it  must  be,  because  Holland's  extended  formula  for 

astronomical  refraction  neglects  the  term  tan3(z).  Compare 

formula  (l6)  of  this  paper. 


However,  the  error  coefficients  in  the  last  two  columns  of 
Table  1  in  [Case  (1962)]  are  insignificant,  as  are  the  similar  values 
in  [Jones  (1961)] .  The  results  given  in  Table  2  in  [Case  (1962)]  , 
labelled  "Flat-Earth  Satellite  Refraction,"  lead  to  correct  values 
for  satellite  refraction,  within  the  limits  imposed  by  the  approxi- 
Hs 

mation  d  =  GQs  ^  ,  namely,  flat  earth  geometry  and  assuming 

Ha  =  0.  The  effect  of  a  flat  earth  geometry  can  be  judged  from 
Table  VI,  columns  (5)  and  (6)  in  this  report. 

The  results  obtained  in  this  paper  indicate  that  proper  consid¬ 
erations  of  refraction  corrections  for  both  the  geodetic  satellite 
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triangulation  and  the  topographic  photogrammetric  satellite  method 
are  possible  without  any  computational  difficulties.  Furthermore, 
it  seems  justified  as  long  as  the  camera  or  the  target  point  is  out¬ 
side  of  the  effective  atmosphere,  to  state  that  refraction  cannot  be 
made  the  scapegoat  if  in  the  future  the  need  should  arise  to  explain 
discrepancies  in  the  results  of  photogrammetric  space  triangulation 
methods . 

This  statement  should  not  detract  from  the  problem  of  scintil¬ 
lation,  which  affects  particularly  short -duration  electronic  flash 
photography.  The  geometrical  significance  of  directions  measured  to 
such  targets  is  limited,  independent  of  the  precision  of  the  measur¬ 
ing  method,  by  the  scintillation  effect  to  an  accuracy  (one  sigma 
level)  of  +£  to  +3  seconds  of  arc. 
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